A channel Brownian pump powered by an unbiased external force 
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Abstract 

A Brownian pump of particles in an asymmetric finite tube is investigated in the presence of 
an unbiased external force. The pumping system is bounded by two particle reservoirs. It is 
found that the particles can be pumped through the tube from a reservoir at low concentration to 
one at the same or higher concentration. There exists an optimized value of temperature (or the 
amplitude of the external force) at which the pumping capacity takes its maximum value. The 
pumping capacity decreases with increasing the radius at the bottleneck of the tube. 
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I. INTRODUCTION 



Noisy transport far from equilibrium play a crucial role in many processes from physical 
and biological to social systems. Molecular motors are the paradigm of how to extract 
useful mechanical energy in the random environment of thermal fluctuations^ - -. Inspired 
on Feynman's ratchet and pawl machine, the appearance of directed motion enhanced by 
thermal fluctuations has been named after the ratchet effect 4 . 

The idea of applying the ratchet mechanism to model pumps has already appeared in 
the literature^— . Prost and co-workers 5 studied the transport of an asymmetric pump with 
a simple two-level model and quantified how vectorial symmetry plus dissipation creates 
a macroscopic motion, even in the absence of any externally applied gradient. Astumian 
and Derenyi- investigated a chemically driven molecular electron pump in which charge can 
be pumped through a tiny gated portal from a reservoir at low electrochemical potential 
to on at the same or higher electrochemical potential by cyclically modulating the portal 
and gate energies. Kosztin and Schulten 7 studied the fluctuation-driven molecular transport 
through an asymmetric potential pump and three transport mechanisms: driven by potential 
gradient, by an external periodic force and by nonequilibrium fluctuations. Nonadiabatic 
electron heat pump was investigated by Rey and coworkers 8 . They presented a mechanism 
for extracting heat metallic conductors based on the energy-selective transmission of elec- 
trons through a spatially asymmetric resonant structure subject to ac driving. Wambaugh 
and co-workers^ studied the transport of fluxons in superconductors by alternating current 
rectification. They found that at nonzero temperatures, the asymmetric geometry of the 
ratchet sawteeth automatically converts applied ac inputs into a net dc motion of fluxons. 
In their another work—, they studied how inertia and particle pair interactions can induce 
transverse rectification of particle flows driven through square arrays or chains of triangu- 
lar defects perpendicularly to their symmetry axis. Recently, Sancho and Gomez-marin^ 1 - 
presented a model for a Brownian pump powered by a flashing ratchet mechanism. The 
pumping device was embedded in a finite region and bounded by particle reservoirs. Their 
emphasis is on finding what concentration gradient the pump can maintain. 

The previous works on pump were on the consideration of the energy barriers. The present 
work is extend to the study of the molecular pump to the case of the entropic barriers. 
The entropic barriers may appear when coarsening the description of a complex system for 



2 



simplifying its dynamics^ - —. We emphasize on finding how the Brownian particles can be 
pumped through an asymmetric tube from a particle reservoir at low concentration to one 
at the same or higher concentration in the presence of an unbiased external force. Our other 
task is on investigation what factors can affect the pumping capacity. 



II. GENERAL ANALYSIS 

We consider overdamped Brownian particles moving in an asymmetric finite tube [Fig. 
1] in the presence of an unbiased external force. The tube is embedded in a finite region and 
bounded by two particle reservoirs. Its overdamped dynamics is described by the following 
Langevin equations^ 1 ^ 

7§ = f(i) +&(*), (1) 

7^ = <.(*), (3) 
where x, y, z are the three-dimensional (3D) coordinates, 7 is the friction coefficient of the 
particle. £ x ,y,z{t) are the uncorrelated Gaussian white noises with zero mean and correlation 
function: < £i(t)£?(£ ) >= 27 kBT5ij5(t — t') for i,j = x,y,z. T is the absolute temperature. 
ks is the Boltzmann constant. < ... > denotes an ensemble average over the distribution 
of noise. Imposing reflecting boundary conditions in the transverse direction ensures the 
confinement of the dynamics within the tube. The shape of the tube is described by its 
radius, 

cj(x) = a sin(— — ) + b, xq < x < xq + L, (4) 

where a is the parameter that controls the flat degree of the tube of the tube, L is the length 
of the tube. The radius at the bottleneck is = b — a. xq is the coordinate of the left end. 
F(t) is an unbiased external force and satisfies 

F(()= K nr< t <nr + Ir; (g) 
[ -F , nr + ^r <t<(n + l)r, 

where r is the period of the unbiased force and F is its magnitude. 
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When the particles move in a confined media, its movement equation can be described by 
the Fick- Jacobs equation^— which is derived from the 3D (or 2D) Smoluchowski equation 
after elimination of y and z coordinates by assuming equilibrium in the orthogonal directions. 
The complicated boundary conditions of the diffusion equation in irregular channels can 
be greatly simplified by introducing an entropic potential that accounts for the reduced 
space accessible for diffusion of the Brownian particle. Reduction of the coordinates may 
involve the appearance of the entropic barriers and an effective diffusion coefficient. When 
\u (x)\ < 1, the effective diffusion coefficient reads^"— 

where Dq = ksF '/7 and a = 1/2 for three dimensions. The prime stands for the derivative 
with respect to the space variable x. 

Consider the effective diffusion coefficient and the entropic barriers, the dynamics of a 
Brownian particle moving along the axis of the 3D tube can be described by the equation^^ 

where a free energy A(x,t) = E — TS = —F(t)x — Tkehihfa) is defined^ - —: here E = 
—F(t)x is the energy, S = fcsln/i(x) is the entropy, h(x) is the dimensionless transverse 
cross section tt[u(x)/L} 2 of the tube in three dimensions. j(x,t) is the probability current 
and p(x, t) is the particle concentration. 

If F(t) changes very slowly with respect to t, namely, its period is longer than any other 
time scale of the system, there exists a quasistatic state. In the steady state, the concentra- 
tion is just a function of space thus the flux becomes a constant j. The concentration p{x) 
follows a first order non homogeneous linear differential equation, whose formal solution is 

, fi'W,,, . f x dz . r z A'(y) , 

Though, unlike typical Brownian motors, the boundary conditions are not periodic nor 
the normalized condition is imposed, the unknown constant cq can be found by imposing 
the left reservoir concentration p Q = p(x ) and the right concentration pi = p(x + L) as 
fixed boundary conditions^. Then c = p(x ) and 

j(Fo ) = Wo - Pl e-ffi] 

LT '[1 + oo'(x)^e~^^[^dx 



The average current is 

J=- f T mt))dt= 1 -[j(F )+j(-F )]. (10) 
t Jo z 

For studying the pumping capacity, we consider the situation in which J tends to zero 
which corresponds the case in which the pump in maintaining the maximum concentration 
difference between the two reservoirs across the membrane with no net leaking of particle. 
The method is the same that in Ref. 11. This situation is analogous the stalling force in 
Brownian motors. 

From Eqs. (8-10), we can obtain 



Fi)(x — xn) Fr\(x — Xr\) . . 

C +L [lW(xffle^^ +e-^-][t*$] 2 dx 

" (11) 



III. RESULTS AND DISCUSSIONS 



For simplicity, we take k B = 1, 7 = 1 and 
is valid only when \uj (x)\ < 1 and its validity 
Thus, we take a = j- and |w'(x)| << 1 which 



L = 2tt throughout the study. Our method 
is testified by using numerical simulations 16 , 
ensure the validity of our method. 



A. Current J at p% = po 

In order to check the effect of the magnitude for the tube, we plot the current J as a 
function of period numbers N( we have N cells in the channel) in Fig. 2. It is found that 
the current J is independent of the period numbers N. Each cell has the same structure 
and pumping capacity. Thus, our results are not finite size effects. 

Figure 3(a) shows the current as a function of the initial coordinate Xq for different values 
of Fq at po = Pi- I n this case, the transport behavior is similar to that in a periodic tube. The 
current is positive for < xo < | and §7r < xo < 2tt, zero at xo = | and §7r, and negative 
for I < xo < |7r. The positions of the negative peak and positive peak are determined 
by the parameters of the system. In order to illustrate the transport behavior, currents 
j(F ), — j(Fo) and J versus xq for F = 0.5 are shown in Fig. 3(b). When < xq < f and 
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§tt < x < 2ir, j(F ) is more than —j(—F ), so the current J is positive. When | < a; < §7r, 
— j(— F ) is more than j(F ) and J is negative. j(F ) is equal to — j(— F ) at x — | and 

§7T. 

In fact, the sign of the current is determined by the asymmetry of the tube. The corre- 
sponding shape of the tube is given by Fig. 4. It is obvious that the initial coordinate xq 
determines the asymmetry of the tube. When the left space of the tube is bigger than the 
right (xq = 0), the particles run to the right in the presence of the unbiased external force. 
The particles move to the left when the right space is bigger (x = n). There is no current 
for the symmetric tube (x = |, |7r). Hence, the system can pump the particles from the 
left to the right for < xq < f and |7r < x < 2n. We take x = in the rest study. 

B. Ratio pi/po at J = 

We now focus on the maximum concentration ratio at both end of the tube. The analytical 
results are given by Eq. (11). In Figs. 5-7, we investigate the pumping capacity and see 
how the concentration ratio depends on temperature T, the radius r b at the bottleneck and 
the amplitude Fq of the external force. 

Figure 5 presents the concentration ratio pi/po as a function of temperature T at x = 0. 
The curve is observed to be bell shaped. When T — > 0, the particles cannot reach all space 
of the tube and the effect of the asymmetric entropic barrier disappears and there is no 
pumping capacity (pi/po = 1)- When T — > oo, the effect of the unbiased external force 
disappears and the pumping capacity goes to zero, also. Therefore, there is an optimized 
value of T at which the ratio takes its maximum value, which indicates that the thermal 
noise can facilitate the particles pumping. 

The ratio pi/po as a function of the radius r b at the bottleneck is plotted in Fig. 6. When 
r b — > 0, no particle can pass through the tube, thus, pi/po — > oo. With increasing of the 
radius r&, the effect of the asymmetric entripic barrier reduces and the pumping capacity 
decreases. 

Figure 7 shows the ratio pi/po vs the amplitude F of the unbiased external force. When 
F Q — > 0, only the effect of the asymmetric entropic barrier exists, so the system has no 
pumping capacity. The ratio pi/po saturates to 1 in the large amplitude F limit. Therefore, 
the pumping capacity is maximum when F takes an optimized value. 
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C. Current J at p\ > po 

In Fig. 8, the current J versus temperature T is presented for different values of p\ at 
p Q — 1. For lower values of pi the current is larger. When pi > p , the current is negative 
for both too low or too high temperature. When T — > 0, the effect of the asymmetric 
entropic barrier disappears, the concentration difference dominates the transport, thus, the 
current is negative. When T — > oo, the effect of the unbiased external force disappears, the 
transport is dominated by the concentration difference and the particles move to the left. If 
the concentration difference exceed its maximum value, the current will be always negative. 

Figure 9 shows the current J as a function of the radius r b at p — 1, Pi — 4. The curve is 
observed to be bell shaped. When r& — > 0, the tube is jammed and the particles cannot pass 
through the tube. When r b — > oo, the effect of the asymmetric entropic barrier disappears, 
and the current is dominated by the concentration difference, thus the current is negative. 
So there exists an optimized value of r b at which the current takes its maximum value. 

IV. CONCLUDING REMARKS 

In this study, we study the transport of Brownian particles moving in an asymmetric 
finite tube in the presence of an unbiased external force. The pumping device is embedded 
in a finite region and bounded by two particle reservoirs. It is found that the particles can be 
pumped from low concentration reservoir to higher concentration reservoir. When p\ = p , 
the current is positive for < xq < ix/2 and 2>n/2 < xq < 2it, zero at xq = n/2 and 37r/2, 
negative for 2ir < xq < 3tt/2. The sign of the current is determined by the asymmetry of 
the tube. The maximum ratio of the concentrations (J = 0) at both ends of the tube is also 
studied. We can find that there is an optimized value of temperature T (or the amplitude 
Fq of the external force) which gives the maximum ratio of pi/po- The ratio pi/po decreases 
with increasing the radius r b at the bottleneck of the tube. We also study the current J as 
a function of J and r b . We observe that the current can tend to its maximum value when T 
or r b take the optimized value. 

Though the model presented does not pretend to be a realistic model for a biological 
pumps, the results we have presented have a wide application in many systems, such as Na, 
K-ATPase pumps^, membrane proteins 18 , and electron pumps^. 
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V. CAPTION LIST 



Fig. 1. Scheme of the pumping device: a spatially asymmetric tube is embedded in 
a finite region of length L and bounded by two particle reservoirs of concentrations po 
and p\. The shape of the tube is determined by its radius u>(x) and the left end coordi- 
nate xo of the tube. The particles in the tube are powered by an unbiased external force F(t). 

Fig.2. Current J as a function of period numbers N at x = 0, a = b = a = 1/2, 
F = 0.5, T = 0.5, po = 1, and Pl = 1. 

Fig. 3. (a) Current J as a function of x for different values of F at a = 6 = 
a = 1/2, T = 0.5, po = 1, and pi = 1. (b)Currents j(F ), — j(— F ) and J versus Xo for 
F = 0.5. The other parameters are the same as that in (a). 

Fig. 4. The tube shape for different values of x . The tube is symmetric at x = n/2 and 
3tt/2. The left space is bigger than the right at x = and the right space is bigger at x = it. 

Fig. 5. Ratio of concentrations as a function of T at a = b = a = 1/2, F — 0.5, 
and xo = 0. 

Fig. 6. Ratio of concentrations as a function of r b at a = ct = 1/2, F = 0.5, T = 0.5, 
and xo = 0. 

Fig. 7. Ratio of concentrations as a function of F at a = 6 = a = 1/2, T = 0.5, 
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and xq = 0. 

Fig. 8. Current J vs temperature T for different values of pi at a = b = a = 1/2, 
F = 0.5, x = 0, and p — 1. 

Fig. 9. Current J as a function of r b at a = a = 1/2, F = 0.5, T = 0.5, rr = 0, 
p = 1, and pi = 4. 
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